Photon-photon correlations and entanglement in doped photonic crystals 
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We consider a photonic crystal (PC) doped with four-level atoms whose intermediate transition is 
coupled near-resonantly with a photonic band-gap edge. We show that two photons, each coupled 
to a different atomic transition in such atoms, can manifest strong phase or amplitude correlations: 
One photon can induce a large phase shift on the other photon or trigger its absorption and thus 
operate as an ultrasensitive nonlinear photon-switch. These features allow the creation of entangled 
two-photon states and have unique advantages over previously considered media: (i) no control 
lasers are needed; (ii) the system parameters can be chosen to cause full two-photon entanglement 
via absorption; (iii) a number of PCs can be combined in a network. 



PACS number(s): 42.50.Gy, 03.67.-a, 42.70.Qs 
I. INTRODUCTION 

Nonlinear effects whereby one light beam influences an- 
other require large numbers of photons |Q or else photon 
confinement in a high-Q cavity p|. Hence the impedi- 
ment towards constructing quantum logical gates oper- 
ating at the few-photon level. The ingenious attempt 
to achieve increased photon-photon coupling in a gas by 
means of control laser fields || has resulted in cooper- 
atively enhanced single-photon absorption and emission 
(as for excitons in solids), but not in two-photon entan- 
glement [Q. A promising avenue has been opened by 
studies of enhanced nonlinear coupling via electromag- 
netically induced transparency (E1T) in gases in the pres- 
ence of control laser fields, which induce coherence be- 
tween atomic levels || . These studies have predicted the 
ability to achieve an appreciable nonlinear phase shift us- 
ing extremely weak optical fields or a two-photon switch 
in the TV-configuration of atomic levels HJ7J . Further im- 
provement of the sensitivity of these schemes has been 
suggested using a rather involved system, in which a sec- 
ond species of coherently driven A-atoms provides the 
matching of the group velocities of interacting photons 

Here we point out that photon-photon nonlinear phase 
shifters and switches are realizable with a potentially 
very high efficiency and without external laser fields in 
photonic crystals (PCs) |9|-[lH| dilutely doped with four- 
level atoms. These atoms have two transitions tuned to 
the two incident photons and an intermediate transition 
tuned to a singular feature of the structured density of 
modes (DOM) spectrum of the PC. The strong nonlin- 
ear effects analyzed in Sees. || and III arise from the co- 
herent Autler-Townes splitting of atomic emission lines 
near a photonic band-gap (PBG) edge |p^-[l5|| and the 
strong photon-atom coupling via photonic defect modes 
in the doped PC [jl6|,[l7). Transparency near a band edge, 
which has been previously predicted for an atomic three- 
level A-configuration is shown here (Sec. Ill) to be 
most suitable for two-photon absorption switching. This 
mechanism, which has not been studied in PCs, is demon- 
strated to be considerably more efficient than its coun- 



terpart in Ref. JtJ. It is predicted to allow complete ab- 
sorption of one photon in the presence of another photon, 
and thereby the creation of a fully entangled two-photon 
state (Sec. [Pvj ). By contrast, appreciab le ( but limited) 



III) to arise by 



nonlinear phase- shifting is shown (Sec. 
tuning the photon frequency to the Raman resonance 
with the defect mode frequency associated with atomic 
doping. These features have unique advantages for quan- 
tum information applications (Sec. IV), as compared to 
previously considered media. 



II. THEORY 
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FIG. 1. (a) Photonic crystal dilutely doped with atoms lo- 
cated at black dots, (b) Four-level atom coupled to a struc- 
tured continuum {A} near the band-edge or defect mode 
frequencies (DOM plotted) via the intermediate transition 
|2) — s- 1 3) and interacting with two photons ui a and u>b at the 
sideband transitions |1) — > |2) and |3) — > |4), respectively. 

We examine the nonlinear coupling of two optical fields 
E a = f a £a& and — fb£bb, where a and b are the re- 
spective annihilation operators, £\ — (hu>i/2eoVi) 1 ^ 2 is 
the field amplitude per photon Ui (i — a, b) within the 
quantization volume Vj = o~ih, and fi(z,t) is the single- 
photon wavepacket envelope. The two fields propagate 
along the 2-axis in a PC dilutely doped with identical 
four-level atoms. The level configuration of the atom and 
the DOM are depicted in Fig. ffl, where the unperturbed 
atomic levels \j), j = 1 . . . 4, have the corresponding en- 
ergies Hujj . The incident photons at frequencies uj a and Ub 
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interact with the atoms via the transitions |1) — > |2) and 
|3) — > |4), respectively, while the transition |2) — > |3) is 
coupled to the structured PC mode-continuum {A}. In 
the dipole and rotating-wave approximations, the Hamil- 
tonian of the Zth atom+ficld can be written as 



contains a 

term 7 4 |Sl h | 2 /(A 2 . + 74) resulting in an additional relax- 
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where fcj (i = a, 6, A) is the wave number of the corre- 
sponding mode, gi — Hjkfi£i/h its atom-field coupling 
strength, fj,jk being the dipole matrix element for the 
atomic transition \j) —> \k), (rv£ = \j){k\ is the atomic 
operator, c A is the A mode annihilation operator and 
J2\ ~ ¥ J doj\p(ui\), where p(uj\) is the DOM of the struc- 
tured continuum. 

We assume that initially the two incident photons 
are in the product state |l a ) 1 1&) , the atoms are in 
the ground state |1) and the continuum is in the vac- 
uum state |0a) - Then the wave function of the system 
reads = A x |1, A , 1„, l b ) + A 2 |2, A , a , U> + 

E A ^3, A |3,1 A O a ,l b )+E A A 4 ,a|4,1a,0 q ,0 6 ) . With the 
Hamiltonian (|l|), the Schrodinger equation leads to the 
following set of equations for the slowly-varying (during 
an optical cycle) probability amplitudes Ay. 



the right-hand side of Eq. (Eq) 

fc i 2 /(Ar ■ — 

ation of the amplitude A-$^\ and thus destroying the co- 
herence between levels |1) and |3) which induces the ab- 
sorption of the uj a photon. To minimize this decoherence, 
we take |At| ^> |f2&|,74, so that the depletion of can 
safely be neglected. In the opposite case of small detun- 
ing | Aft | < 74, it has been suggested to use this effect for 
constructing a sensitive two-photon switch |jj . 

Under these conditions, one obtains effectively free- 
space propagation of the E b field. By contrast, the Rabi 
frequency of the E a field, in the slowly varying envelope 
approximation, obeys the following propagation equation 
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il a = iaVl a , 
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with the solution fl a (z,t) = fi a (0, i — z/v g )exp(iaz). 
Here the macroscopic complex polarizability a is given, 
under the weak-field linear-response assumption, by 
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(2a)is the linear resonant absorption coefficient on the atomic 
transition |1) — > |2), with oo the resonant absorption 

(2b) cross_sect ^ on ' N the density of doping atoms and n a the 
(averaged) refraction index for the Lu a photon, and 
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where tti = p,jk{0i \ E$ \ li)/h = gi is the Rabi frequency of 
the corresponding field at the position of the Zth atom, 
A a = uj a - u 2 i, A b — uj b - uj i3 and A A = lox ~ w 23 
are the detunings from the respective atomic transition 
frequencies, and jj (j — 2,3,4) is the relaxation rate of 
level \j), which accounts for both spontaneous radiative 
decay at a frequency far from PBG so that it is treated as 
a Markovian process (Wigner-Weisskopf approximation) 
and nonradiative (mainly vibrational) relaxation in the 
PC. 

The Wigner-Weisskopf approximation is inapplicable 
near the PBG where the DOM varies rapidly [Q. Thus, 
in the vicinity of the transition frequency u> 23 , we must 
exactly integrate the last term of Eq. ( pb| ) for the spe- 
cific PBG model employed. To solve Eqs. (g), we make 
the weak-field approximation A\ ~ 1 (much less than one 
photon per atom) and use the second-order perturba- 
tion theory thus obtaining the steady-state expressions 
for the atomic response. In doing so, one can see that 
the probability amplitude A^^x is inversely proportional 
to the detuning A b . Then, near the Raman resonance 



is the integral of the saturation factor over the structured 
DOM. The group velocity v g is expressed as 
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At the exit from the medium z = the delay time T^ e \ of 
the field, relative to the passage time Tq — (,n a /c through 
a passive medium, can be written as 
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del 
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The slowly varying field-propagation equation (|^) must 
be justified upon examining the group- velocity dispersion 



D = 



d 2 k a d 2 Re(a) 
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which is responsible for the spreading and reshaping of 
the photon pulse. In analogy with the weakly interacting 
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Bose gas, D can be interpreted as being inversely propor- 
tional to the photon "mass" . In what follows, conditions 
such that D is small are discussed for ui a . Due to the 
weak- field approximation A\ ~ 1 and large detuning 
the u)b photon propagation is nearly free and its group 
velocity dispersion is negligible. 

To calculate integral (||), we assume the isotropic PBG 
model |l2|-|l5f with the atoms doped at the positions of 
the local defects in the PC separated by a distance d from 
each other. These "impurities" of the crystal structure 
form the defect modes in the PBG, which are localized at 
each atomic site in a volume Vd — (rL) 3 of several (r) 3 
lattice cells L 3 Q^Q. In the dilute regime d > rL, 
one can neglect dipolc-dipolc interactions and tunneling 
("hopping") of photons between the atoms Jb|[l{| and 
the defect modes can serve as high-Q cavities. For W23 
near (or within) the PBG frequency, we have L ~ ttc/lu23- 
Hence, the dilute regime limits the dopant density to 
N < (a>23 / 'ncr) 3 . Then, in the vicinity of the upper edge 
uj p of the PBG, the DOM function can be written as 



p(uj) = PdS(uJ - LOd) + Pi 



Q(lj - to e ) 



(9) 



where Q(lu) is the Heaviside step function, pd and p e are 
PC-specific constants ]l^ , ^5| , and u>d is the frequency of 
the defect mode. The integration of Eq. (g), with p(u>\) 
given by Eq. (H) and A& 3> A 0i x, leads to 
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(10) 

where Ad. e = uJd.e — ^23 <C ^23 are the detunings of the 
defect-mode and PBG-edge frequencies from the atomic 
73 +ji\n b \ 2 /A 2 b the |1) «-> |3) 



resonance UI23, 731 = 
decoherence rate, S3 
level 1 3), and 



= |fib| 2 /Af, the ac Stark shift of 
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are the coupling constants of the atom with the struc- 
tured reservoir, whose main contributions are near uid 
and u». 



we plot the spectrum for two different values of the cou- 
pling constants (3d and (3 e , that there is, however, a sub- 
stantial difference between the spectra in the foregoing 
frequency regions, for the following physical reasons: 
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FIG. 2. (a),(b) Imaginary (solid lines) and real (dashed 
lines) part of the complex polarizability a, as a function of 
the detuning A a for the case fit = 0, ao = 1 cm" 1 , = — 1, 
A e = 1 and 73i = 0.001. (a) /3 d = p e = 1; (b) (3 d = (3 e = 2. 
(c) Delay time Tdei (P er un it z), and (d) group velocity dis- 
persion coefficient D as a function of the detuning A a for the 
same parameters as in (a). All parameters are in units of 72. 

o) In the vicinity of A^, the radiation emitted by the 
atom at the frequency uid remains confined in the de- 
fect mode for a long time as in a high-Q cavity. As 
the u> a photon wave packet approaches the Ith atom, 
^a//?d( < C 1) fraction of its amplitude is transferred into 
the defect mode, inducing the corresponding population 
of level 1 3) (assuming adiabatic Raman transfer |2C[|). At 
the end of the wavepacket, when Sl a — s- 0, all the popu- 
lation of level 1 3) returns to |1) and the radiation that 
has been confined to the defect mode is added to the 
tail of the propagating wave packet, until it encounters 
the next atom. This strong interaction of the atom with 
the defect mode at A a ~ A^ splits the spectrum by the 
amount equal roughly to 2 fid and causes EIT ||. The 
transparency window is rather broad and is given by the 
inverse Lorentzian [see the first term on the right-hand 
side of Eq. (tHl)]. The corresponding group velocity is 
much smaller than the speed of light: 



III. RESULTS 

To illustrate the results of the foregoing analysis, we 
first plot in Fig. || the polarizability (Q), delay time (0) 
and group- velocity dispersion (g) for the case of one inci- 
dent photon uj a (P-b = 0). Clearly, two frequency regions, 
A a ~ Ad and A a ~ A e , where the absorption vanishes 
and, at the same time, the dispersion slope is steep, are 
of particular interest. One can see in Fig. 0(a), (b), where 



<9Re(a) 



duj a 



72 ao 



< c , 



(11) 



which leads to a large delay time T c i c i = Cl v g [Fig- ||( c )]- 
One has to keep in mind, however, that the absorption- 
free propagation time is limited by the EIT decoherence 
time Tdoi < 73] 1 j5j, which imposes a limitation on the 
length C of the active PC medium. The correspond- 
ing group- velocity dispersion (^) is small (a large photon 
mass), and, therefore, will not cause much spread of the 
u a wave packet [Fig. 0(d)]. 
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b) In the vicinity of A e , the strong interaction of the 
atom with the continuum near the band edge u> e causes 
the Autler-Townes splitting of level |2) into a doublet 
with a separation equal roughly to (3 e . One compo- 
nent of this doublet is shifted out of the PBG while 
the other one remains within the gap and forms the 
so called photon-atom bound state ]T^- [T5| , p^ | . Conse- 
quently, there is vanishing absorption and rapid variation 
of the dispersion at A a ~ A e . Here the delay time can 
be large, but the group velocity dispersion is also very 
large [Fig. |2|(c),(d)]. Hence, as uj a is tuned very close to 
the band edge, the transmitted pulse shape is distorted. 
Since the transparency region is very narrow with a width 
Aoj ~ 73 (<§; j2,0e), f° r an absorption- free propagation 
of the Lo a photon, the temporal width of its wavepacket 
i~ a = l a /c should satisfy the condition r a > tt/Alo. Simul- 
taneously, a small deviation from the condition A a = A e 
will lead to a strong increase in the absorption of the E a 
field. 



1.2 



' 

(a) 


i 




'\\ i 1 i 1 i 










A A°' 3 


\ 1 1 




/ l i 7 \ o 
'/ 7 \ C 

/ M '/ \ 

/ \ 1 \ 
/ l\ 1 \\ 






0.5 1 1.5 2~ 


i , i \\y i ^ 






33 



1 

0.8 
0.6 
0.4 
0.2 

0.9 

0.6 

0.3 



-0.3 
-0.6 



FIG. 3. (a) Imaginary, and (b) real part of the complex 
polarizability a as a function of the detuning A a for the case 
fit, — (solid lines) and fib = gb (dashed lines). All pa- 
rameters are the same as in Fig. ^(a) and S3 = —0.1 (i.e., 
At — — 10|fi(,| 2 ). The insets magnify the important frequency 
regions. 

Let us now switch on the Uh photon. As seen from 
Eq. (|l0|), its effect is merely to shift the spectrum by the 
amount equal to S3 (Fig. [|). This shift, however, will 
have different implications in the two frequency regions 
distinguished above: 

i) If S3 -C (3d, he., the Stark shift is smaller than the 
width of the transparency window at A^, the medium 
will still remain transparent for a Lu a photon with the 
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detuning A a = Ad, but its phase will experience an 
preciable nonlinear shift <p a given by 



dRe(a) 

(pa = Rc(a)z ~ s 3 z 

OLU a 



72 a 



1% 



(12) 



ii) On the other hand, for a u a photon with the de- 
tuning A a = A e , the medium, which is transparent for 
fib = 0, having 



j / x 72^0 13 ^ 



(13) 



will become highly absorptive (opaque) even for such a 
small frequency shift as S3 (provided | S3 1 > Aw), chang- 
ing (13) to 



Im(a) 




(14) 



and thus acting as an ultrasensitive, effective switch. 
These are the main results of the present work. 

The remaining question is how to maximize the inter- 
action between the uj a photon, which propagates with 
a small group velocity (|ll[), and the u>b photon, which 
propagates with a velocity close to the speed of light. A 
possible technique to achieve equal group velocities for 
both photons, as suggested in Ref. @], is to have a sec- 
ond kind of A-atoms in the interaction region and apply a 
driving held that would result in EIT and reduced group 
velocity for the tjj\> photon. To keep our scheme simple, 
we have chosen here not to adopt a similar approach. 

In our system, the interaction between the photons 
is maximized if: they enter the medium simultaneously; 
the transverse shapes of their tightly focused (<r a ,fc ~ o~o) 
wavepackets overlap completely; and the quantization 
(wavepacket) length lb of the ujb photon satisfies the con- 
dition (lb + ()/c < Q/vg. Then the u>t, photon leaves 
the medium not later than the u> a photon. The effective 
interaction length between the two photons is therefore 
z cff ~ h v g/ c 5~ Cj after which the two wavepackets slip 
apart. This effective interaction length has the following 
implications: 

1) Since s 3 cx 1/4, the interaction- induced phase shift 
( |l2l ) saturates with distance and is the same for all If, 
satisfying the above condition: 



dRe(a) lbV g 
duj a c 



|M34.| 2 ^b 

2eo?iC(TfcA(, 



(15) 



2) As was mentioned in Sec. || and carefully analyzed 
in Ref. 0, in the case of a small detuning |A(,| < 74, the 
presence of the u>b photon will result in the destruction 
of EIT in the vicinity of A Q = A^ and, consequently, the 
strong absorption, which is determined by 
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With the effective interaction length z c g ~ lf,v g /c, where 
v g is given by Eq. (|ll]), the power loss at the exit from 
the PC is 



2Im(a)z eff : 



(16) 



Thus, in the vicinity of the two-photon (Raman) res- 
onance with the defect mode, A Q = Ad, the absorp- 
tion ( |l6| ) (| A;, | < 74), as well as the phase shift jl5| ) 
(|A;,| 3> , "74,), saturate over a distance equal to z e g. 

3) The dramatic advantage of the present scheme over 
conventional EIT schemes ||,|7j is that, in the vicinity of 
the two-photon resonance with the band edge frequency, 
A a = A e , the group velocity of the uj a photon wavepacket 
is close to the speed of light in the presence of the nearly- 
free propagating uib photon wavepacket [see Fig. 0(c)]. 
Therefore, there is no velocity mismatch of the two pho- 
tons and no saturation of the absorption with distance, 
the interaction length being close to the length £ of the 
PC. Thus the absorption probability of the ui a photon, 
1 — exp[— 2Im(a)£], where Im(a) is given by (|l4|), can be 
made arbitrary close to unity by choosing a long enough 
PC. 

It can be checked that, for the parameter values used 
in Fig. |, and ^34 ~ 1 a.u., w 6 = 4x 10 15 rad/s (A& = 470 
nm), |A 6 | = 1072 , 72 ~ 5 x 10 7 s~\ <r b ~ lO" 10 cm" 2 , 
we obtain the power loss 2Im(a)£ ~ 0.46£ at A a = A e 
and the phase shift <fi a ~ 0.1 rad at A a = A<j. Thus, the 
presence of one u>b photon induces either strong absorp- 
tion or a large phase shift of the u> a photon, depending 
on the frequency region employed. 



wherein the states in which the uo a photon is or is not 
absorbed are equally superposed. 

b) In the case of the ui a photon being initially in a state 
with a rather well defined phase (i.e., coherent state) \a) 
and having the detuning A a ~ Ad, after passing through 
the PC, the resulting state of the system is given by 

-L(|10) W6 + 101),.*)® \a) Ua 

-> i( \10) ut \e*+*a) Ua + \01U \a) Ua ) , (18) 

wherein the states in which the tu a photon does or does 
not acquire the phase shift 4> a are equally superposed. 
The state ( |l8| ) is fully entangled only if \e l ^ a a) and \a) 
do not overlap in the phase plane (the same requirement 
as for a "Schrodinger cat" state in a PC |l7|]). 

These features can be used to appreciably advance to- 
wards the goal of producing entangled states of radiation 
or logical photon switches for quantum information pro- 
cessing, owing to the unique advantages of the doped 
PCs over conventional EIT schemes ^-g] or high-Q cav- 
ities 0): (i) No control lasers are needed to create the 
quantum interference responsible for the EIT effects, (ii) 
The system parameters can easily be adjusted to provide 
full two-photon entanglement via absorption, (iii) Sev- 
eral doped PCs can be combined via dispersive couplers 
either to enhance the accumulated nonlinear phase-shift 
or to perform a chain of logical operations. The fact that 
the band structure of dilutely doped PCs is insensitive to 
the locations or exact concentration of the dopants and 
defects should facilitate the incorporation of several PCs 
in one network. 



IV. DISCUSSION 

The foregoing results have demonstrated the ultrahigh 
sensitivity of photonic absorption or phase shift in a 
doped PC to the presence of an additional photon. Using 
the procedure of Ref. one can employ these features to 
construct two types of entanglement between the photons 
Lo a and uib' 

a) Suppose that a photon pair \l) bJa \l)u b is simulta- 
neously generated in some parametric or other process. 
The ui a photon enters the doped PC. The u>b photon is 
split between the two arms of a 50%-50% beam splitter, 
resulting in an entangled state 1/V2(|10) W( , + |01) Wb ). 
One arm of this state, say 1 10)^ , enters the PC together 
with the ui a photon, while the other arm, |01) Wij , does 
not. In the case of the ui a photon with the detuning 
A a ~ A e , after passing through the PC, the resulting 
state of the system has the fully entangled form 



|01)c 



II) 



-^(|io) w jo) Wa + |0i) w ji)O 



(17) 
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